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ABSTRACT 


Previous measurements of parametric x-radiation (PXR) spectra from mosaic 
graphite observed yields that greatly exceeded theoretical predictions. At the Bragg 
scattering condition, the interaction of relativistically charged particles and the atomic 
planes of a crystal produces scattering of virtual photons. This is known as PXR. The 
discrepancies between actual yields and theoretical predictions might possibly have been 
caused by the mosaic spread of the atomic planes of the crystal. Re-analysis of previous 
data has shown a miscalculation of the effective target thickness used in earlier work. 
There was also an error in the calculation of the theory that was used for comparison. 
Correction of the two errors brought the data and theory in closer agreement, but the 
mosaic spread does change the intensity of the PXR. 

In this thesis, three samples of thick mosaic carbon graphite crystal with differing 
mosaic spreads are compared to determine the effect of the mosaicity on PXR yield. By 
obtaining “rocking curves” where the target orientation angle is changed with respect to 
the incident electron beam and observation angle, it was found that the mosaic spread does 
have an effect on the PXR spectral distribution. This effect is quite significant and can be 
seen when viewing the “rocking curves” of the samples. As the mosaicity increased, the 
angular distribution was spread out and the coherence condition was relaxed. 




vi 




TABLE OF CONTENTS 


I. INTRODUCTION 1 

II. THEORETICAL BACKGROUND 3 

m PXR EXPERIMENTAL DESCRIPTION 9 

A. ELECTRON BEAM ALIGNMENT AND STEERING 9 

B SPECTROSCOPY ELECTRONICS 11 

C. ENERGY CALIBRATION.12 

IV DATA AND ANALYSIS 19 

A. CALCULATION OF PXR PEAK AREAS 19 

B CALCULATION OF PXR ABSOLUTE YIELD 28 

C. CALCULATION OF INTENSITY RATIOS 32 

D. CORRECTIONS TO PREVIOUS DATA 42 

E ERROR ANALYSIS 48 

V. CONCLUSIONS AND RECOMMENDATIONS 57 

APPENDIX A SAMPLE THEORETICAL YIELD PROGRAM 59 

APPENDIX B. SAMPLE ERROR ANALYSIS PROGRAM 65 

APPENDIX C DATA VERSUS THEORY FOR ALL SAMPLES 67 

LIST OF REFERENCES 75 

INITIAL DISTRIBUTION LIST 77 


vii 





















viii 




ACKNOWLEDGEMENTS 


First and foremost, we would like to thank our thesis advisor X. K Maruyama, 
for his patience which we tried on a daily basis. Without his expertise and guidance, this 
work would not have been possible. We would also like to remind him that great 
innovations come from stupid mistakes, and we made many of those. 

We would also like to thank Don Snyder and Harold Rietdyk for the long, late 
hours they spent with us during the data collection process We also need to thank Don 
for repairing all of the damage we did to the computer and software during our work, 
including the worn out mouse. 

Q. Li of the Naval Surface Warfare Center provided us with a base Mathcad 
program to work from for which we are very grateful Without this, the theoretical 
calculations would have been very painful. We also need to express our thanks to Don 
Rule, also of Naval Surface Warfare Center for all of the information and calculations that 
he provided to us. Also, Mel Piestrup of Adelphi Technology sacrificed of himself and his 
time to help us and we greatly appreciate that We would also like to wish him good luck 
with the new data he is collecting. 

Finally, we would like to thank our families: Karen. Christopher. Alyssa, Aaron, 
Max, and Trey for putting up with the long hours and our absence in the past year. For 
this, we dedicate this work to them. 


IX 




I. INTRODUCTION 

Scattering of virtual photons produced by the interaction of a relativistic charged 
particle and a crystalline material at or near the Bragg condition, 2dsin0=nX, results in the 
production of parametric x-radiation (PXR). [Ref. 1] This scattering is associated with 
the Coulomb field of the relativistic particle and the atomic planes of the crystal; in this 
case, a thick mosaic carbon graphite crystal was used. When the incident electron beam 
angle with respect to the crystal axis satisfies the Bragg condition, virtual photons are 
placed on the mass shell and scattered as real photons. The first PXR theory was 
developed by Ter-Mikaelian (1971) as resonant radiation that was produced in a thin 
crystal. [Ref. 1] Development of the theory for a thick crystal followed with experiments 
in the mid 1980’s. [Ref. 2] In the early 1990’s the first PXR experiments outside of the 
former Soviet Union were performed at Naval Postgraduate School (NPS), Monterey, 
California. The first experiments using a carbon target at NPS were reported in June 1992 
[Ref. 3] and in December 1992. [Ref. 4] These experiments used a compression annealed 
pyrolytic graphite (CAPG) with a mosaic spread of 0.45°. 

The recent theoretical work has been done by Rule et al at the Naval Surface 
Warfare Center in Silver Spring, Maryland. [Ref. 5] This theory predicts experimental 
yields due to the influence of the mosaic spread in an attempt to explain the discrepancies 
seen between previous experiments and theories. 

In this thesis, thick mosaic graphite crystals with mosaicities of 0.45°, 1.31°, and 
2.5° were used with a 90 MeV electron beam at a Bragg angle of 22.5° which yielded x- 
rays in the 4 to 30 keV range. Analysis of these spectra confirmed that mosaic spread 
does indeed have an effect on the angular distribution and coherence conditions. 
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H. THEORETICAL BACKGROUND 

PXR is an important source of x-rays because of the high energy photons that are 
produced with the relatively low electron energy requirements of the experiments. PXR is 
produced by the interaction of relativistically charged particles and the atomic planes of a 
crystal. When the Bragg condition, 2dsin0=nX, is satisfied, virtual photons are placed on 
the mass shell and scatter out as real photons. Ter-Mikaelian was the first to produce a 
theory for PXR in a thin crystal. The first theory for a thick crystal was by Feranchuck 
and Ivanshin. [Ref. 6] The emitted photons spectral and angular distributions can be 
expressed by: 


— =—oiL\ I 2 
dn ± dco 2jt a 


K r ,COV + t 


’K-*x) 2 +£(l-v 2 )] 2 


-exp(--£-)]%). 


0) 


where <7 = ^-- yjo) 2 -K 2 + x z —y [Re(g 0 ) - 6 S ] (2) 

where N is the number of photons per electron, e is the electron charge, 0 S is the multiple 
scattering angle, L a =[coIm(go)]’\ go and g t are the perpendicular and parallel Fourier 
components of the dielectric susceptibility, x is the reciprocal lattice vector, n x is the 
photon index of refraction, to is frequency of the emitted photon, v is particle velocity, k x 
is the component of the wave vector perpendicular to the crystal plane. Feranchuck and 
Ivanshin also found that the number of photons per electron per solid angle per electron 
energy can be found using angular displacement as follows: 

d^N_ 1 e 2 1 1 , 2 e x 2 cos 2 (2e g ) + e/ _ 1 _ 

de,3e y 3co 2x'hca, g '' + 

where 0 X and 0 y represents the angular displacement from the Bragg condition, where p s is 
defined as X/[L a (co)27t] with L a (co)is the photon absorption length, (Ob is the frequency of 
the emitted photon at Bragg condition , y is defined as l/Vl-(v/c) 2 , Aco is the difference 
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between the emitted photon © and © B . [Ref. 7] To obtain the angular distribution of the 
spectral PXR spectra. Equation (3) can be integrated over ©. 

This theory when used by DiNova [Ref. 4] yielded values of N that were much 
smaller than the experimental photon production rate for all spectral orders and for each 
angle at or near the Bragg condition for the “near field” case and values that were larger 
than experimental yields for the far field case. Table 1 shows a sampling of values from 
DiNova and the theoretical values produced by the above equation integrated over ©. 
There were corrections made to the values calculated by DiNova in the value of the 
effective tin thickness as will be presented in the data analysis section. 



Near Field 

Far Field 

Angle (0 B ) 

Data 



Data 

EHSI 

Theory w/ 
Scattering 

20.1° 

3.80e-6 

2.53e-6 

2.49e-6 

- 

- 

- 

21.1° 

1.36e-5 

5.47e-6 

5.25e-6 

2.27e-6 

2.059e-6 

1.99e-6 

21.6° 

3.01e-5 

1.70e-5 

2.21e-5 

5.46e-7 

1.15e-5 

9.03e-6 

22.1° 

4.70e-5 

2.18e-5 

1.51e-5 

1.14e-6 

1.18e-5 

7.85e-6 

22.6° 

6.22e-5 

2.18e-5 

1.38e-5 

4.68e-6 

8.57e-6 

4.949e-6 

23.1° 

3.82e-5 

2.18e-5 

1.5le-5 

3.87e-6 

1.18e-5 

7.85e-6 

23.6° 

1.58e-5 

1.70e-5 

2.21e-5 

8.74e-7 

1.15e-5 

9.03e-6 

24.1° 

8.07e-6 

5.47e-6 

5.25e-6 

4.25e-7 

2.059e-6 

1.99e-6 

24.6° 

6.47e-6 

4.26e-6 

4. : . e-6 

3.06e-7 

1.23e-6 

1.21e-6 

25.6° 

4.76e-6 

2.53e-6 

2.49e-6 

- 

- 

- 


Table 1. Comparison of Past Experimental Data Corrected 

with Theoretical Calculations. 


'or Effective Tin Thickness 


Rule’s calculations presented at Tomsk by Rule et al [Ref. 5] convolves PXR 
spectral-angular distributions with Gaussian models of the mosaic tilt angle distributions 
and the effects of mosaic spread, a m , is brought into play. The effect of multiple 
scattering, a s , is also included in the following: 
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and the effects of mosaic spread, c m , is brought into play. The effect of multiple 
scattering, cr s , is also included in the following: 



with g x 2 = < 7 y 2 = <j m 2 + a s 2 used to add the mean square mosaic tilt angle to the mean 
square scattering angle. 0 xl , 6^, 0 yl and 0 y2 are the limits of integration and in this case 
are all equal to the radius of the circular detector aperture in radians, (Oe(n) = mtc/dsin(0 B ) 
where d is the interplanar spacing, y is the Lorentz factor, % 0 is the mean dielectric 
susceptibility, and Xio is the dielectric structure factor This theory appears to more 
accurately support the experimental results obtained. 

The effects of finite aperture size must be considered when working with PXR. 
These effects depend on the aperture area and the distance between the target and the 
detector. From the above, the solid angle of view for the detector can be determined and 




this can be used to determine how much of the PXR is actually being viewed by the 
detector. 

There are two regions for a finite detector aperture called the narrow line limit 
regime and the aperture limited bandwidth regime. [Ref. 4] For the first case, A6 x =ps,. 
For the second case, A0 x »p s and the distribution can be approximated as 

— - Ae,y,(«) for -e7<e,<e7 (! 

N 

and is equal to zero otherwise. Where 



and a y - - 2L where v ' 

V ph 

0 p h = m 2 /E 2 + 0 S 2 +1 go I ■ This is the regime with which this thesis will be concerned. 

Within this regime, the “near field” is defined as A0 X > FWHM of J 2 (u) and the 
“far field” is defined by A0 X < FWHM of J 2 (u). The “near field” and “far field” cases are 
presented in Figures (1) and (2), respectively. For this thesis, the experiments were done 
in the “near field”, however, there are some references to previous “far field” experiments 
and results. As can be seen in Figure (1), in the “near field” there is very little finite 
aperture effects and most of the PXR is within the field of view of the detector. On the 
other hand, care must be taken in the analysis of the “far field” data as a large portion of 
the PXR is out of the view of the detector due to its limited solid angle size. 
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Figure 1 . h and the field of view for the “near field” 



Figure 2 . J 2 and the field of view for the “far field” case 
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HI. PXR EXPERIMENTAL DESCRIPTION 

The Naval Postgraduate School Electron Linear Accelerator (LINAC) was used to 
accelerate dark current electrons to energies of approximately 95 MeV. [Ref. 7] These 
electrons then struck three separate targets of compression annealed pyrolytic graphite 
(CAPG) of differing mosaicities. Initial energy calibration was performed using 
sandwiched foils of tin (Sn) and Yttrium (Y). In order to provide in situ energy calibration 
and the ability to calculate photon yields [Ref 4], a tin (Sn) foil was placed behind the 
samples. “Rocking curve” data was obtained by placing the target so that its crystal plane 
was placed at an angle of 21° with respect to the electron beam direction and then rotating 
the target in 0.5° increments from the electron beam until 26°. This allows for comparison 
of the relative intensities of the n* order peaks between samples. 

A. ELECTRON BEAM ALIGNMENT AND STEERING 

The LINAC is shown in Figure 3. All operations are performed and monitored 
from the control room. Remote cameras are used in the end station to assist in obtaining 
proper electron beam alignment. The accelerator consists of three 10 ft. sections, each 
with a separate Klystron connected by waveguides. The LINAC is operated at 2856 MHz 
with a puise repetition rate of 60 Hz. The beam macro structure length is approximately 1 
M-sec. Thus, the LINAC duty cycle requires the ability to count at 60 MHz when a single 
photon is detected during each LINAC pulse. [Ref. 4] In order to satisfy the requirement 
for low count rate, the LINAC was adjusted during operation to maintain one count every 
five to ten machine macropulses. This also ensured that simultaneous counts of two 
photons would not appear as one count of a higher energy photon. This is of concern as 
PXR peaks are integral multiples of each other with respect to energy and could result in 
“pile up.” 

A reference laser and retro optics is used to position the target ladder on which all 
target samples are mounted. A beam is then established in the accelerator after a vacuum 
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is drawn in the target chamber. The quadrupole magnets are used to provide coarse 
steering of the electron beam towards the target Fine adjustment of beam position is 



Figure 3. Overhead view of the Naval Postgraduate School 100 MeV Linear Accelerator 
(LINAC). 

obtained by steering the beam through a small pinhole in the target ladder using the end 
station steering magnets 

Target chamber geometry is shown in Figure 4. The alignment laser is located 90° 
from the entrance port. The Canberra Si(Li) detector [Ref. 8] is located 135° from the 



Figure 4. Overhead view of LINAC end station experimental setup. 
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entrance port. The target ladder could be raised and lowered remotely and it consisted of 
a two axis goniometer to allow rotation of the mosaic crystal. This same detector 
geometry was maintained for all data runs. 

B. SPECTROSCOPY ELECTRONICS 

As with all detectors, the Canberra Si(Li) detector is very susceptible to ground 
loops as well as radiated noise from the klystrons. Thus, the end station is enclosed in a 
metal mesh cage to reduce the interference of radiated noise from other sources. Lead 
and borated polyethylene bricks were also used around the Si(Li) detector to reduce the 
background. Figure 5 shows the electronics setup. All system components were located in 
the control room with the exception of the detector and amplifier which were located in 
the LINAC end station. 



TENNELEC 

SifLf) detector h ^|TC 244 
Amplifier 


TENNELEC TC 
■/308 Dual Linear 
Gate 



Figure 5. Electronics setup. 

The signal received by the Si(Li) detector is amplified by a Tennelec TC 244 amplifier. A 
Tennelec TC 308 dual linear gate was used to reduce unwanted signal by providing a 
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25 nsec gate during which signal could be received. This effectively made the counting 
system only active during the LIN AC machine pulse. Triggering of the gate was 
performed using a Stanford Research System Inc. model DG 535 digital delay/pulse 
generator. By introducing a delay time (T D ) with respect to the machine s start sequence 
time (T 0 ). the gate was adjusted to coincide with the arrival of the beam pulse, thus, T D 
was adjusted to 46 psec. 

Another problem inherent in energy calibration of the spectra is the affect of 
klystron noise on the signal. The LINAC produces 60 pulses of one microsecond duration 
every 60 seconds. This gating has no effect on the calibration spectrum and assists in 
limiting noise which is non-coincident with PXR and Sn fluorescence production. On the 
other hand, klystron noise adds a negative bias to the received detector signal which 
produces a negative bias to the calibration curve if one relies solely on radioactive sources 
for energy calibration. Figure 6 shows the klystron noise signal in relation to detection 
system gating For this reason calibration spectra were taken prior to each day's data run 
using Sn and Y to eliminate this potential source of error 



Figure 6. Representative Shape of Klystron Noise During PXR Data Collection. Line 

Indicates Approximate Gated Portion. 


C. ENERGY CALIBRATION 

During the initial energy calibration conducted each day prior to the rocking 
curves the LINAC was operated at full gun grid voltage. The beam strikes sandwiched 
foils of tin (Sn) and Yttrium (Y) in order to observe the and Kp fluorescent lines. 
Energy of the Y Ka is 14 957 keV and the energies of the Sn Ka and Kp lines are 
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25.270 keV and 28.483 keV respectively. [Ref. 9] PCA-II pulse height analyzing 
software [Ref. 10] was used to ensure calibration of all 512 channels. “Rocking curve” 
data was then taken beginning at 21.0° from the electron beam and proceeding in 0.5° 
increments to 26.0° with the LINAC operating at dark current. This is required in order 
to limit the photon count rate since the photon detector, a Canberra model SI200250 
Si(Li) detector, has a resolution time of 12 psec determined by the pulse shaping 
preamplifier. Dark current is the term used to describe operation of the LINAC with no 
gun grid voltage applied. Thus, only stray electrons are accelerated toward the target 
crystal. 

In order to provide in situ energy calibration, a tin (Sn) foil was placed behind each 
CAPG sample. Table 2 provides target and tin (Sn) thickness for each sample. Each 
thickness measurement was obtained by taking the average of five thickness 
measurements at different positions on the sample. Each target crystal has a density (p) of 
2.260 ± 0.005 grams per cubic centimeter. The {002} reflecting plane spacing is 3.356 ± 
0.003 Angstroms. 

The Union Carbide specifications for the mosaicities of the samples were 0.4±0.1°, 
1.2±0.2°, and 3.5+1.5°. Since the tolerances from the prescribed mosaic spread were 
quite large, the samples were sent to Naval Surface Warfare Center (NSWC) for 
determination of the actual mosaic spreads. This was measured with an x-ray 
diffractometer. The resulting spectra from the x-ray diffraction are presented in Figures 7 
through 9. [Ref. 11] For each sample, the mosaic spread was determined by the frill width 
at half maximum (in degrees) of these curves. The resulting spreads were 0.45°, 1.2°,and 
2.5°. 


Mosaic Spread 

0.45° 

! 

1.31° 

2.5° 

Target Thickness 

1389 ± 1 pm 

1659 ± 1 pm 

1789 ± 1 pm 

Sn Thickness 

152 + 1 pm 

154 ± 1 pm 

156 ± 1 pm 


Table 2. Sample and Sn Thickness for each Sample. 
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_ Energies (in keV) for each channel were established during the calibration runs by 
the PCA-II pulse height analyzing software using the following equation: 

E = ao + ai(channel) + a 2 (channel) 2 (12) 

where ao, a.\ and a 2 are constants determined by performing Gaussian fits of the counts 
received in each region of interest surrounding the Yt and Sn calibration peaks. These 
constants for each series of data is shown in Table 3. The energy calibration was verified 
by comparing the K a and K g lines produced by the Sn backing used during each PXR data 
run against the known value. Good agreement was observed for all PXR spectral data. 


Mosaic Spread 

0.45° 

1.31° 

2.5° 

a 0 

1.95 

1.86 

1.34 

ai 

0.0722 

0.0712 

0.0764 

a 2 

0 

0 

8.53e-6 


Table 3. Calibration Constants for each Sample. 
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Figure 8. Mosaicity Measurement curve for 1.31° Sample. 


16 











r 


18 



IV. DATA AND ANALYSIS 


A. CALCULATION OF PXR PEAK AREAS 

The data from the detector was collected using a PCA-II pulse height analyzer 
[Ref 10] in ASCII format. Raw data spectra for the 0.45°. 1.31°. and 2.5° mosaic 
samples are shown in Figures 10. 11. and 12 respectively. Data was recorded for each 
sample spectrum for approximately 30 minutes with electron energies (E e ) of 99 05 MeV, 
95.32 MeV. and 97.62 MeV respectively. The spectra were then input into Quattro Pro 
(QPRO) spreadsheets [Ref. 12] so that it would be in proper format for the spectrum 
analyzing software. Peakfit v3.0 [Ref 13] was used to remove the background noise and 
produce a Gaussian fit for each spectral peak. Figure 13 a) b) and c) show the 
uncorrected, background removed, and attenuation corrected spectrums respectively for 
the 0 45° sample at 22.5°. The program also integrated each peak to produce a net area 
for each peak and the full width at half maximum values for each peak of the spectrum. 
The values were placed in a Microsoft Excel worksheet [Ref 14] and corrected for 
attenuation using Photcoef. [Ref 15] Photcoef was used to calculate the attenuation 
coefficients for the Kapton (C22O5N2H10) window, the Berylium (Be) window and the air 
gap at each energy. Thicknesses for each of these were 0 0025 cm 0.005 cm. and 0.5 cm 
respectively. At low energies, the attenuation coefficients varied across each peak. For 
the first peaks only, the attenuation coefficients for the left most channel, the peak center 
and the right most channel were added in quadrature and averaged to obtain a more 
accurate correction for the peak. For the higher order peaks, the attenuation coefficients 
were approximately constant across the peak and the coefficient was obtained for the 
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center channel of the peak and applied to the entire peak. The areas were then corrected 
to account for these attenuations as follows: 

^ _ jy * gl(f l P t )B'! + (MP t )A.r + (MP t )Kapion] ( 13 ) 

where No is the total number of photons before attenuation. N is the total number of 
photons after attenuation, p. is the attenuation coefficient, p is the density of the 
attenuating medium, and t is the distance travelled in the attenuating medium. The 
efficiency of the detector had to also be accounted for in the corrections The efficiency 
curve for the Canberra detector at various energies can be found in Figure 14. [Ref. 8] 
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Figure 11. Raw 1.31° Mosaicity Spectra at E e = 95.32 MeV 
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Figure 13. a) Uncorrected PXR spectrum for 0.45° mosaic spread at 22.5°. b) The same 
spectrum with background removed, c) Spectrum corrected for attenuation between the 
sample and the detector. 
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Figure 14. Efficiency Curves for Canberra Si(Li) Detector. 

B. CALCULATION OF PXR ABSOLUTE YIELD 

The corrected PXR areas are then used as follows to calculate the number of 


photons produced for each electron. 




(Area PXR / \ 

^ Atime J 


Combining Equations 14 and 15 yields: 


Area PXR cr 

TIN 

47T_ 


(Area m )AW n 
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where AreapxR and Area-nN are the integrated intensities of the respective peaks, corrected 
for both attenuation and detector efficiency, AWtin is the atomic weight of tin c is the 
electron interaction cross section (1 Oe-26 m 2 ) [Ref 16] Q det is the detector solid angle 
(4.15e-4 sr), f dex is the de-excitation transition probability (0.712) [Ref. 17], and t is the 
effective tin thickness. The effective thickness used in the above equation was calculated 
by dividing the thickness of the tin (ti) by the sine of the angle of incidence of the electron 
beam. The thickness of the tin (ti) was 152 + 1 pm, 154 ± 1 pm, and 156 ± 1 pm for the 
0.45°, 131°, and 2.5° mosaic samples respectively. The geometry of the electron beam as 
it strikes the Sn foil presented in Figure 15. This calculation was done for the first peak at 
each angle of the rocking curve for all three samples The number of photons per electron 
(N) versus angle for the first and second order PXR peaks are presented in Figures 16 and 
17 for each sample. 



Figure 15. Configuration in Target Chamber Used to Determine Tin Effective Thickness, 
ti is the Normal Target Thickness. 0 B is the Bragg Angle or the Angle Between the 
Target Face and the Electron Beam Direction. 
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Figure 16. Comparison of Corrected First Order PXR Intensities (N). 0.45°. 1.31°, and 2.5° Refer to the Mosaic Spread of the 
Graphite Targets. The Electron Energies (E e ) are 99.05 MeV. 95.32 MeV, and 97.62 MeV Respectively. 


N(photons/electron)x 10 A -7 
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SECOND PEAK COMPARISON 
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C. CALCULATION OF INTENSITY RATIOS 


The higher order peaks were compared to the first peak absolute yield for each 
angle of the samples to give an intensity ratio. These calculations are presented in Table 
4 and will be used for analysis of theoretical predictions. For each sample, the values 
presented in Figure 16 were normalized to the highest first order intensity to provide a 
comparison of the relative intensities of the first order PXR peaks across the “rocking 
curve.” These values are presented in Figure 18. The same comparison was made for the 


second order PXR peaks and is presented in Figure 19. 

In order to provide a comparison of the data with theory for all three samples, 
Equation 3 was integrated over ©, 0 X , and 0 y for the solid angle subtended by the detector. 
These calculations are referred to here as theory without scattering. To account for 
multiple scattering and mosaic spread. Equation 17 was used. [Ref. 4] These calculations 


d~ N _ e~ Q) l 

se,de,~h 4 *“ 




where S tt 2 =p-+8/+|* 0 | 

2Ild 0, “ = 0^muilscai 0^beam-div 4 m 0^n 


sin* d B Q’ +9 y ' + Q pk ~ 


(17) 

(18) 
(19) 


are referred to as theory with scattering. Appendix A shows a sample of the MathCad 
[Ref. 18] program used to calculate the values shown for theoretical calculations, with and 
without taking into account the effects of mosaic spread and multiple scattering. Figures 
20 through 25 give a comparison of the data with the resulting theoretical curves for each 
sample, including both the first and second order PXR peaks. Agreement between the 
data and theory appears to break down as the mosaic spread increases. Thus, the method 
presented by Rule et al of convolving the PXR spectral-angular distributions with 
Gaussian models of the mosaic tilt angle appears to provide better agreement. Table 5 
provides the theoretical intensities and ratios for the first and second order PXR peaks at 
Bragg condition provided by Rule for a detector solid angle of 4.15e-4 sr. [Ref. 19] 
Further calculations using this model are in progress. 
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FIRST PEAK NORMALIZED 





SECOND PEAK NORMALIZED 
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Figure 20. Corrected First Order PXR Intensity Comparison with Theory for 0.45° Mosaic Sample. 
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0.45 DEG. FIRST PEAK DATA VS, THEORY 
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Figure 21. Corrected Second Order PXR Intensity Comparison with Theory for 0.45° Mosaic Sample. 























Figure 22. Corrected First Order PXR Intensity Comparison with Theory for 1.31° Mosaic Sample. 
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1,31 DEG, FIRST PEAK DATA VS. THEORY 




1.31 DEG. SECOND PEAK DATA VS. THEORY 
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Figure 23. Corrected Second Order PXR Intensity Comparison with Theory for 1.31° Mosaic Sample. 




Figure 24. Corrected First Order PXR Intensity Comparison with Theory for 2.5° Mosaic Sample. 
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2.5 DEG. FIRST PEAK DATA VERSUS THEORY 
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Figure 25. Corrected Second Order PXR Intensity Comparison with Theory for 2.5° Mosaic Sample. 







Mosaic Spread 

N = 1 

N = 2 


0.45° 

2.14e-5 

7.81e-6 

0.36 

1.31° 

1.63e-5 

5.79e-6 

0.36 | 

2.5° 

1.06e-5 

3.70e-6 

0.35 


Table 5. Theoretical Intensity Comparisons Using Rule et al for Detector Solid Angle of 

4.15e-4 sr. [Ref. 19] 

D. CORRECTIONS TO PREVIOUS DATA 


The same 0.45° mosaic spread sample was used in recent experiments as was used 
in prior experiments. DiNova [Ref 4] measured photon yield (N) with the detector placed 
at 100 cm and 29 cm from the target which will be referred to as the “far field” and “near 
field” cases respectively. DiNova calculates the effective thickness of tin used in the 
absolute yield calculations by dividing the thickness of the tin sample by the cosine of 0b- 
[Ref. 16] As presented earlier, this should be divided by the sine of 0 B . Once this 
correction is made to DiNova’s data, the comparison between her data, and recent data 
reveals similar results. Differing detector solid angles account for most of the difference 
between recent data and the older near and far field data. Other differences include 
varying tin foil thicknesses, different detectors, and different data analysis procedures. 
Table 6 summarizes these differences. In order to determine how large the effect of the 
differing data analysis techniques was, DiNova’s raw data was reprocessed in the exact 
manner as for the three samples presented here. For comparison purposes, Table 7 shows 
the first order peak data using the corrected effective tin thickness from DiNova’s thesis 
and DiNova’s data re-analyzed. Figures 26 and 27 compare this data graphically. An 
attempt was made to re-analyze the data for the 25.6° spectrum, but the first order peak 
was not distinguishable. DiNova’s re-analyzed near field and far field and the recent 0.45° 
mosaic sample data are also con ired in Table 8 to show how the distance from detector 
to target effects the absolute yield. It is important to note that the new data fell at a 
distance that was in between DiNova’s far and near field results. 
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Near Field Data 

New Data 

Far Field Data 


2.39e-3 sr. 

4.15e-4 sr. 

2.0e-4 sr. 

Tin Thickness (tj.) 

| 



Detector 



■EEIOSIESH 

Analysis Procedure 

Summed Counts in 
PXR Peak 

Peakfit v3.0 

Summed Counts in 
PXR Peak 


Table 6. Summary of Differences between Prior and Recent Experiments. 


Angle (0b) 

Near Field 
Data 

Near Field Data 
Re-Analyzed 

Far Field 
Data 

HM 

20.1° 

3.80e-6 

3.60e-6 

- 

- 

21.1° 

1.36e-5 

1.28e-5 

2.27e-6 

5.12e-7 

21.6° 

3.01e-5 

2.73e-5 

5.46e-7 

9.30e-7 

22.1° 

4.70e-5 

2.71e-5 

1.14e-6 

2.55e-6 

22.6° 

6.22e-5 

5.49e-5 

4.68e-6 

8.49e-6 

23.1° 

3.82e-5 

3.30e-5 

3.87e-6 

2.07e-6 { 

23.6° 

1.58e-5 

1.28e-5 

8.74e-7 

1.59e-6 

24.1° 

8.07e-6 

7.33e-6 

4.25e-7 

7.75e-7 

24.6° 

6.47e-6 

4.77e-6 

3.06e-7 

4.59e-7 

25.6° 

4.7 6e-6 

- 

- 

- 

Table 7. Comparison of DiNova’s 

Data Corrected for t 

re Effective Tin Thickness and 


DiNova : s Data Re-Analyzed Using PeakFit (0.45° Mosaic Spread) 


By changing the distance from the target to the detector, the solid angle of the 
detector is changed. As can be seen in Figure 28 for the three different solid angles, there 
is a marked effect on the photon yield (N) as the distance is increased 
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1 

Near Field Data 

New Data 

Far Field Data 


2.39e-3 sr. 

4.15e-4 sr. 

2.0e-4 sr. 


N(photons/elec.) 

N(photons/elec.) 

N(photons/elec.) 


3.60e-6 



l 

|N> 

b 

0 


1.21e-6 


-1.5° 

1.28e-5 

1.50e-6 

5.12e-7 

1 

® 

0 

2.73e-5 

3.29e-6 

9.30e-7 

-0.5° 

2.71e-5 

7.81e-6 

2.55e-6 

0° 

5.49e-5 

2.27e-5 

8.49e-6 

0.5° 1 

3.29e-5 

1.63e-5 

2.07e-6 

1.0° 

1.28e-5 

9.12e-6 

1.5 6e-6 

1.5° 

7.33e-6 

3.99e-6 

7.75e-7 

2.0° 

4.76e-6 

2.85e-6 

4.59e-7 

2.5° 


1.56e-6 


O 

® 


1.18e-6 



Table 8. Comparison of 0.45° Mosaic Sample Results by Solid Angle. 


As can be easily seen from Figures 16 and 17 as the mosaic spread of the crystal 
increases, the angular distribution of the “rocking curve” also increases. The PeakFit 
software was again used to obtain the full width at half maximum for the first order peaks. 
The resulting distribution appears linear and is shown in Figure 29 Thus, increasing 
crystal mosaicity relaxes the coherence condition and spreads out the resulting rocking 
curve” angular distribution. 
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Figure 27. Corrected and Re-Analyzed Data for Far Field 0.45° Mosaic Sample for Detector Solid Angle of 2 Oe-4 sr. 
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Figure 29. Full Width at Half Maximum of First Order Peaks as a Function of Mosaic 
Spread. 

E ERROR ANALYSIS 

There are two possible sources of errors in this experiment, namely, systematic and 
random errors. The sources of systematic errors are caused by limitations of the systems 
used and effect all measurements in the same fashion. These errors are often hard to 
analyze and to even detect by statistical analysis Random errors are experimental 
uncertainties that can be revealed by doing the same measurements a number of times to 
get an average value of the measurement. Treatment of the two types of errors are 
handled differently. [Ref 20] 

In handling systematic errors, the best analysis is to predict the possible sources 
and verify that these errors will be less than the precision that is required for the 
experiment. The systematic errors for this experiment are the angle of target incidence 
which is due to the mechanical motor that turns the target and is subject to over and/or 
under rotation. The beam alignment is also a systematic error that would not have an 
effect on each individual run, but when comparing different days runs, this error will be 
evident. This error also effects the number of electrons that actually strike the target 
crystal and tin foil which is a value that is calculated from the areas of the spectra. The 
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combination of these errors produced a shifting of maximum intensities of 0.5° for the 
0.45° mosaic sample, and 1.0° for the 1.31° and 2.5° samples. This error is evident when 
comparing the maximum experimental intensities to the calculated theoretical peak 
intensities. 

The random errors are handled by adding in quadrature the errors involved with 
the variables used in calculating the number of photons per electron as seen in Equation 
16. The only two terms which contribute to this calculation are the areas of the PXR 
peaks and the area of the tin peak. Equation 20 was used to calculate the fractional error 
in the number of photons per electron. The random errors in area were analyzed using 
Equation 21. 



where F is the counts in a spectral peak, ao is the height of the PXR peak, ai is the center 
channel of the PXR peak, a 2 is the standard deviation of the curve fitted peak used by 
Peakfit to determine the best fit peak, and dao, dai, and da 2 are the errors in each of these 
values respectively that are provided by the Peakfit Peak Analysis Program. The integral 
was processed using a MathCad [Ref. 18] program that is listed in Appendix B for each 
peak of the PXR spectrum and the tin peak for each spectrum. The errors in the PXR 
peak and the tin peak are then added in quadrature to obtain the total error for each 
calculation of the number of photons per electron. Table 9 lists the results of Equations 
20 and 21 as well as those from adding in quadrature the errors in the Area P xR and Area™ 
using Equation 22 for each of the first peaks for each spectra. Figures 30 through 34 
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1 

0.45° Mosaicity 

1.31° Mosaicity 

2.5° Mosaicity 

EH5B0S1 

MCAD 

1/VArea 

MCAD 

1/VArea 

MCAD 

1/VArea 

21.0° 

.0079 

.062 

.0078 

.054 

.0164 

.042 

21.5° 

.0082 

.064 

.0091 

.036 

.0192 

.042 

22.0° 

.0113 

.041 

.0107 

.019 

.0247 

.027 

22.5° 

.0262 

.026 

.0252 

.012 

.0515 

.028 

23.0° 

.1747 

.018 

.0499 

.010 

.0405 

.023 

23.5° 

.0805 

.018 

.0782 

.014 

.0411 

.023 

24.0° 

.0303 

.023 

.0420 

.015 

.0345 

.022 

24.5° 

.0131 

.035 

.0131 

.027 

.0276 

.027 i 

25.0° 

.0094 

.035 

.0174 

.036 

.0214 

.033 

25.5° 

.0070 

.052 

.0089 

.057 

.0149 

.037 

26.0° 

.0061 

.059 

.0069 

.062 

.0154 

.039 


Table 9. Estimation of the Random Error in PXR Intensity (N) Calculation Using 


Equations 20, 21, and 22. 

show the random error based on Equations 20 and 21 in the earlier near and far field re¬ 
analyzed 0.45° mosaic spread data as well as the new data for each of the first order 
peaks. Using this method, larger intensity peaks incurred a larger fractional error due to 
the relatively larger size of the peak. 
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Figure 31. Random Error Distribution for 0.45° Mosaic Sample With Detector Solid Angle of 2.0&-4 sr. 
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Figure 34. Random Error Distribution for 2. 5° Mosaicity Sample With Detector Solid Angle of 4 1 !>e-4 sr. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

As suspected from previous experimental and theoretical comparisons, the mosaic 
spread of the crystal used does have a definite effect on the outcome of the PXR spectra 
and the resulting intensities and absolute yields. Much work is being done by Rule et al to 
model theoretically the effect of the mosaic spread on PXR scattering. The theory 
presented by Rule et al [Ref. 5] seems to predict within experimental errors the data 
obtained through the use of the three samples for this thesis more than the older theories 
that were previously used. The work of Rule et al incorporates the mosaic spread, which 
the previous work by Feranchuck and Ivansin [Ref. 6] neglected. As can be seen by 
Figures 16, 17, and 29 an increase in the mosaic spread relaxes the coherence condition 
and spreads out the angular distribution of the scattered photons. Previous data used a 
miscalculated effective thickness for the tin which made the absolute calculated yield 
appear higher than it should have been. In the same comparison, the theory was also 
miscalculated making it to appear lower. By correcting for these two mistakes, the data 
and the theory are more comparable, but there still needs to be the correction for the 
mosaic spread as it does have an effect Theoretical calculations need to be done to 
support the off Bragg condition measurements as well as the higher order peak yields from 
theory presented by Rule et al. Through more experimental data and more theoretical 
calculations to compare with, a better understanding of the mosaic effects on PXR spectra 
can be achieved. 

Future thesis topics should investigate the same samples in a “far field” 
configuration as well as new samples with different mosaic spreads. Contacts should be 
maintained with Rule et al to continue to compare past and present data with new 
theoretical calculations as they become available. 
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APPENDIX A. SAMPLE THEORETICAL YIELD PROGRAM 


Differential cross section calculation 
Carbon (002) 

Input parameters: 

degree 
Angstrom 
gm molar mass 


e B := 22.5 


d := 6.7078 


m a := 12.011 


V := 35.189 


Z a - 6 


cubic angstrom 
Atomic number 


carbon3 

.mcd) 

a :=0.45 

degree 

:= 1659 


a : = 226 

gm/cc 

e ;= 95.32 

MeV 


mos 


1.31*71 

180 


Constants 


e c := 1.6021917* 10 


19 


Electron charge 


m := 9.109558*10' 

c := 299792458 


-31 


a := 7.29735 M0' 


-3 


Electron mass 

speed of light 

Fine structure constant 


6.626196-10 


34 


1 bar " 


(2-tc) 


Plank's constant 


N A •= 6.022169*10 


23 


Avogadro constant 


r e := 2.817938-10' 


15 


Classical electron radius 
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Unit conversion 


e 


e B := K “ 


B 


L = 


L*10 


180 

-6 


sm(0 B ) 


e a = it- 


v = v-io 


180 


30 


Pa :=p a -10 6 d:=d-10- 10 


Bragg frequency calculation 

i .= 4 Order of Bragg reflection 


1‘IZ‘C 


CO 


B 


d-sin(0 g 


Bragg frequency function 


Absorption length calculation 
Absptn := READPRN( carbon) 
i := 0.. rows (Absptn) - 1 
Absptn i 0 := log (Absptn^ 0 ) + 6 


m, 


L a (0»:= 


N A 'Po 


•10 


-linterp(Absptn <0> ,Absptn <1> ,log-— 

e C 


Input is associated with absorption 
data file *.prn generated by xcom.exe 
in units of MeV vs barns 

Absptn^ j := log ('Absptn; ^ - 28 

/(o-ht 


Atomic form factor function 

cf. International Tables for x-ray < 

3] := 2.31 

:« 1.02 % ■- 1-5886 

bj := 20.8439 

b 2 = 10.2075 b 3 := 

c 0 := 0.2156 

i := 1.. 4 


a 4 = 0.865 


f (x) :=^a;-exp| 
i 


-b--lx*10' 


10 


+ c, 


b 4 := 51.6512 


Input parameters 

x should be sin(0)/A, 
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Dielectric function 


r.'(2'it-c) 

r^) =--— 

ic-Voa 


sm (e B ) 

x(co) :=---- CO 

2-n-c 

S £• := 16 Input parameter: Selection factor 


X ioO) = r 5j(©) f(x(a)))- 



% 0 (®) := }( 28 - 8162 )' 



•f(x(©)) 


N calculation 
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56 L = 0.011522 50 A = 0.03-tan(e B ) 

V 50 A 

50 2 =O.O11522 05H 

= 


i = 1..8 


0:=O 

•0 + 50 1 

•50 2 

Z 1 := 


N b ( 0 §* 0 v) d0 § d0 V 


0 - 50 t 

-50 2 


•0 + 50 1 

*50 2 

0 := .008727 ^ := 


N b( 9 ^ 9 ¥ ) d9 ^ de v 


0 - 50 ! 

-50 2 


* 0 + 50 j 

*50 2 

0 := .017453 Z 3 := 


N b( 9 !;’ 0 M/) d 9 i; d % 


0-502 

56 2 


*0 + 50 2 

*50 2 

0 := .02618 z 4 := 


N b( 9 ^ 9 V ) d9 ^ d9 V 


0-50i . 

-56 2 


r 0 + 50, 

r50. 

9 := .034907 

1 

L, 

z 5 = 


N b( 9 £’ 0 y) d 9 S d 9 v|t 


0 - 50 1 

-50 2 


*0 + 50 1 

•50 2 

0 := .043633 Z 6 := 


N b ( 9 £> 9 y) 5® % d9 V 


0-50i . 

-80 2 


•0 + 501 

•50 2 

0:= .05236 


N b( e 5* 9 v) d9 5 d0 V 


0-501 . 

-S 0 2 


•0 + 50 1 

•50 2 

0 := .061087 Z 8 := 




0 - 50 1 

-50 2 
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0.157-Z a -(Z a + 1 2 L 1.13-10 4 (Z a )' 3 '-.511 2 -L 

6msqr 1 ^ . 


>. in n V “ “B (|Z IO^b )!) 1 L a(“ 

N bi( e 5-%) ■——m 5 • 

4-srn e B 


e^ 2 -(cos(2-0 B )) 2 + e v 2 


0.511\ „ 2 . 2 „ /. \2 / 

—— + 0^ + 0^ + ® msqr + (® mos ) + |xo(®B 


i = 1.. 8 


0 = .008727 


0 := .017453 


0 := .02618 


•0 + 50J fS0 2 

0 - 50 j ^-50 2 
1*0 + 50 t I*50 2 


0-60 1 J - 50 2 
'0 + 50 1 1*50 2 




Nu 0 


, ^0 ¥ ) d9 ^ 


N b l 9 


£»%) <*§<*11* 


0-50 1 J-50 2 

'0 + 50 j 1*50 2 


N bi e 4 ,0 v ) d0 ^d0 v 


0 = .034907 


0 := .043633 


0-50J J-50 2 

'0 + 50J rse 2 


N bl( e §. 0 ¥ )*5®y 


0-50J J-50 2 
•0 + 50 2 T50 2 


10-50 2 J-50 2 


N b l 0 


^>0 v ) v 
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0 := .05236 


zz, 


0 = .061087 

zz, 


i:=4 


z = 


z is w/o 


re +■ 58 1 


r 5e 2 


J 0 - 59 j 
r e -h 5e j 


J-502 

r 50 2 


e- 5e : 


-50 2 


N bl ( 0 xj/) d0 ^ d % 

N bl( 0 ^ 0 V ) d0 ^ d0 V 


0 


mos 


1.31-tc 

180 


80 j = 0.011522 
E e := 95.32 
L := 1659 


zz = 


scattering, zz is with scattering 
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APPENDIX B. SAMPLE ERROR ANALYSIS PROGRAM 


inputs for peak height: F is the area under the PXR peak 

a Q := 349.6233 da 0 := 3.7878 F := 1500.9662 

inputs for peak center channel: 
aj := 37.46832 da r = .021196 

inputs for peak width: 
a-, := 1.712699 := .02188 
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APPENDIX C. DATA VERSUS THEORY FOR ALL SAMPLES 

This appendix provides a complete comparison of the actual and theoretical 
intensity ratios between PXR peaks. Theoretical calculations without scattering are based 
on the integration of Equation 3. Theoretical calculations with scattering account for 
multiple scattering and mosaic spread based on the integration of Equation 17 Here we 
also provide a comparison of the observed intensities versus the theoretical intensities both 
with and without accounting for multiple scattering and mosaic spread as calculated using 
the MathCad program found in Appendix A. 
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